Accepted for publication in ApJ 



Imaging Molecular Gas in the Luminous Merger NGC 3256 : 
Detection of High- Velocity Gas and Twin Gas Peaks in the 

Double Nucleus 
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ABSTRACT 

Molecular gas in the merging starburst galaxy NGC 3256 has been imaged 
with the Submillimeter Array at a resolution of 1" x 2" (170 x 340 pc at 35 Mpc). 
This is the first interferometric imaging of molecular gas in the most luminous 
galaxy within 2=0.01. There is a large disk of molecular gas (r > 3 kpc) in the 
center of the merger with a strong gas concentration toward the double nucleus. 
The gas disk having a mass of ~3 x 10 9 M & in the central 3 kpc rotates around a 
point between the two nuclei that are 850 pc apart on the sky. The molecular gas 
is warm and turbulent and shows spatial variation of the intensity ratio between 
CO isotopomers. High- velocity molecular gas is discovered at the galactic center. 
Its velocity in our line of sight is up to 420 km s -1 offset from the systemic velocity 
of the galaxy; the terminal velocity is twice as large as that due to the rotation 
of the main gas disk. The high-velocity gas is most likely due to a molecular 
outflow from the gas disk, entrained by the starburst-driven superwind in the 
galaxy. The molecular outflow is estimated to have a rate of ~10 M yr _1 and 
to play a significant role in the dispersal or depletion of molecular gas from the 
galactic center. A compact gas concentration and steep velocity gradient are also 
found around each of the twin nuclei. They are suggestive of a small gas disk 
rotating around each nucleus. If these are indeed mini-disks, their dynamical 
masses are ~10 9 M & within a radius of 170 pc. 
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1. Introduction 

The luminous infrared galaxy NGC 3256 has an infrared luminosity of 10 1L56 L & at the 
distance of 35.4 Mpc (1"=170 pc; Sanders et al. 2003), which makes it the most luminous 
galaxy within ,2=0.01 (Sargent et al. 1989). The vast luminosity, largely emitted in the far- 
IR, comes from within the central < 20" of the galaxy (Smith & Harvey 1996). The center of 
the merging galaxy hosts a starburst seen across a wide range of wavelengths (e.g., Graham 
et al. 1984; Doyon et al. 1994). High-resolution imaging with the Hubble Space Telescope 
(HST) revealed hundreds of bright young clusters in the galactic center (Zepf et al. 1999; 
Alonso-Herrero et al. 2002). 

The extreme starburst is expected to blow interstellar gas out of the system and may 
leave a gas-depleted elliptical galaxy (Graham et al. 1984). Indeed, an outflow of interstellar 
medium, or superwind, has been detected. Its evidence includes LINER-like optical line 
ratios as well as large line widths (up to 6000 km s _1 ) off the nucleus (Moran et al. 1999), 
blueshifted absorption and emission lines (Heckman et al. 2000; Lipari et al. 2000, 2004), 
and the pattern of optical polarization (i.e., reflection by entrained dust) around the galaxy 
(Scarrott et al. 1996). These observations suggest that the galaxy wind extends out to several 
kiloparsecs. 

NGC 3256 is in the late stage of galaxy merging (Vorontsov-Velyaminov 1959; Toomre 
1977). There is a pair of long tidal tails seen in the optical (see Fig. 1) and in H I emission, 
extending 40 kpc on each side (de Vaucouleurs & de Vaucouleurs 1961; English et al. 2003). 
Such morphology suggests a prograde-prograde merger of two gas-rich spiral galaxies of 
similar size (Toomre & Toomre 1972; White 1979). It also suggests that we view the system 
with a low inclination angle (Feast et al. 1978; English et al. 2003). The radial light profile 
of the galaxy in the f^-band is close to but has not achieved the de Vaucouleurs profile (i.e., 
I(r) oc exp(— At 1 / 4 ) ) of elliptical galaxies (Moorwood & Oliva 1994; Rothberg & Joseph 
2004), unlike the majority of merger remnants surveyed by Rothberg & Joseph (2004). This 
suggests that the merging of the stellar systems in NGC 3256 is not yet complete. 

Likely corresponding to the incomplete merger, two nuclei with ~5" (850 pc on the sky) 
separation have been detected in the near infrared, radio, and X-rays in the starbursting 
center of the merger (Moorwood & Oliva 1994; Norris & Forbes 1995; Lira et al. 2002). 
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They are aligned in the north-south direction, and the southern nucleus is highly obscured, 
rendering it invisible in the optical. Neff et al. (2003) suggested a low luminosity AGN in 
each nucleus from comparison of radio and X-ray observations, while Jenkins et al. (2004) 
found no strong evidence for an AGN in X-ray data alone. 

There is abundant molecular gas (~10 10 M Q ) in the galactic center around the double 
nucleus, probably feeding the starburst (Sargent et al. 1989; Casoli et al. 1991; Aalto et 
al. 1991a; Garay et al. 1993). NGC 3256 is one of the first galaxies in which unusually 
large 12 CO/ 13 CO intensity ratio characteristics of luminous mergers was observed; it was 
attributed to the merger and starburst environment (Aalto et al. 1991b; Casoli et al. 1992b). 

The observations so far are broadly in agreement with the previous observations of 
other luminous and ultraluminous mergers and with the galaxy evolution models involving 
merging. They suggest that a galaxy collision and subsequent merger bring most of gas in 
the progenitor galaxies to the merger center, cause a burst (s) of star formation as well as 
galaxy wind driven by the starburst, and eventually make an elliptical galaxy, sometimes 
with a phase of luminous nuclear activity (Schweizer 1986; Sanders & Mirabel 1996; Genzel 
& Cesarsky 2000; Sanders 2004; Veilleux et al. 2005, and references threin). 

The proximity of a source of this luminosity and the near face-on configuration make 
NGC 3256 an ideal target to test and refine the scenario of merger-induced starburst and 
galaxy evolution. Moreover, NGC 3256 is arguably the southern-sky counterpart of Arp 
220, the archetype of the infrared-luminous merging galaxies (Soifer et al. 1984). The two 
galaxies not only share the title of the most luminous galaxy within their respective distances 
but also the close binary nuclei, since Arp 220, which is twice as distant as and four times 
more luminous than NGC 3256, has two nuclei only 0.3 kpc apart on the sky. Thus they 
are most likely near the final stage of galaxy merger when the merger- induced evolution is 
expected to be most rapid and prominent. Only the southern location of NGC 3256 (Dec. 
= —44°) has delayed one of the most important observations required to better understand 
the merger-starburst evolution — high-resolution observations of its molecular gas. 

We have made high-resolution observations of the cold interstellar medium in NGC 
3256 in 1.3 mm using the Submillimeter Array (SMA) 1 . Our observations aim to uncover 
the distribution, kinematics, and physical properties of the cold molecular gas and dust in 
the merger-starburst environment of NGC 3256 with a high resolution and sensitivity. Of 
particular interest are the distribution and dynamics of gas around the starbursting double 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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nucleus and how the molecular gas is affected by the starburst. The new data of ~1" 
resolution are from the first interferometric observations of molecular gas in this galaxy. 
They provide more than 10-fold improvement in spatial resolution over previous single-dish 
observations made with > 20" beams. 

We begin by introducing our SMA observations (§2), then portray the overall distri- 
bution and kinematics of molecular gas in the central ~9 kpc (§3), report the detection of 
gas concentrations associated with the two nuclei and gas kinematics suggestive of rotation 
around each (§4), report the discovery of high velocity gas and model it as a molecular out- 
flow from the starburst (§5), and characterize the properties of molecular gas in the galactic 
center on the basis of intensity ratios of CO isotopomers (§6). Our observations are discussed 
in the context of the galaxy merger and merger- induced starburst in §7, and summarized in 
§8. 

In this paper, we refer to the north and south nucleus as N and S, respectively, and 
adopt the following positions for them on the basis of radio and X-ray observations (Neff et 
al. 2003; Norris & Forbes 1995; Lira et al. 2002); a(N) = 10 h 27 m 5r!23, 5(N) = -43°54'14"0 
and a(S) = 10 h 27 m 5r!22, S(S) = -43°54'19'. / 2 in J2000. The systemic velocity of the galaxy 
in the literature is V sys (LSR, radio) m 2775 km s _1 (Casoli et al. 1991; Aalto et al. 1991a; 
Garay et al. 1993; English et al. 2003). 

2. SMA Observations and Data reduction 

We observed NGC 3256 at 1.3 mm using the Submillimeter Array (SMA), which con- 
sists of eight 6 m-diameter antennas at the summit of Mauna Kea, Hawaii (Ho et al. 2004). 
We chose the center position of our observations to be a = 10 h 27 m 51!22, 5 = -43°54'19"2 
(J2000), which is centered between the two nuclei. The primary beam of the SMA antennas 
has the FWHM size of 52" (= 9 kpc) at 230 GHz. The extended and compact array config- 
urations were used for our observations in February 20th and April 4th 2004, respectively, 
with all eight antennas participating on both nights. They provided projected baselines 
ranging from 6 m to 179 m. The total integration time on the galaxy was 6.9 hours in two 
tracks. 

The receivers were tuned to simultaneously observe the three CO lines, 12 CO(2-l) in 
the upper side band (USB) and 13 CO(2-l) and C 18 0(2-l) in the lower sideband (LSB). 
The center frequency of the USB was 227.720 GHz, and the LSB was 10 GHz lower. Each 
sideband had 2 GHz of bandwidth and a spectral resolution of 0.81 MHz. The weather was 
excellent on both nights. The 225 GHz zenith opacity measured at the neighboring Caltech 
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Sub millimeter Observatory was 0.05 and 0.08 for the first and second nights, respectively. 
The double-sideband system temperature toward the galaxy was 160 K in median, even 
though the galaxy was observed between the elevation of 14° and 26°. 

We observed the quasar J 1037— 295 once every 20 min for gain calibration. The flux 
density of the quasar was estimated to be 0.86 Jy and 0.84 Jy for LSB and USB, respec- 
tively, from comparison with Mars, for which we assumed the brightness temperature of 
200 K. The system passband was calibrated by observing Jupiter, Saturn, and a few bright 
quasars. The elevation-dependent attenuation by the atmosphere was corrected by using 
the system temperature. The precision of our gain calibration was assessed by comparing, 
between the two tracks, the total 12 CO flux detected in the same range of hour angle on the 
common baselines in the two array configurations. The flux agreed within 10%. We therefore 
assign, conservatively, ±10% for the random error of our flux measurements reported in this 
paper. The error in our absolute flux scaling, however, can be larger due to unknown or 
uncharacterized systemic effects. 

The SMA data were reduced using MIR, which is an IDL version of MMA (Scoville 
et al. 1993), MIRIAD (Sault et al. 1995), and the NRAO AIPS package (Bridle & Greisen 
1994). After the standard passband and gain calibration, the spectral channels that did not 
have emission lines were combined to make continuum data of 0.33 GHz bandwidth in each 
sideband. Channels containing the high-velocity emission that we discovered were excluded 
from the continuum data. The continuum was subtracted from spectral line data in the uv 
plane. Images of the line and continuum emission were then made by Fourier transforming 
the respective set of visibilities. Various spatial and spectral resolutions were achieved by 
changing -uf-weighting and channel binning. Among our datasets are the lowest-resolution 
one made with natural weighting and 20 km s -1 resolution to detect faint extended emission, 
and the highest resolution 12 CO maps made with the super-uniform weighting and 10 km s -1 
resolution. Intermediate resolution maps with the robust uv weighting are also used. The 
rms noise in the natural weighting data cubes is 20 ± 2 mJy beam -1 . The SMA maps 
presented in this paper are not corrected for the primary beam attenuation, except for Fig 4 
(d). The attenuation is approximately Gaussian and is a factor of 2 for a source 26" from the 
map center. The correction for it is made in all flux and spectrum measurements, though it 
is rather small for most emission that we detect. Velocities in this paper are measured with 
respect to the local standard of rest (LSR) and are defined in the ratio convention. 

We compared our data cubes with single-dish observations in the literature, and found 
that our interferometric observations recovered almost all of the flux in the region. The obser- 
vations at the 15 m Swedish-ESO-Submillimeter-Telescope (SEST) gave the total 12 CO(2-l) 
flux of (3.2 - 3.6) xlO 3 Jy km s -1 and 13 CO(2-l) flux of (1.3 - 1.5) xlO 2 Jy km s -1 in their 
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beam at the galactic center (Casoli et al. 1991; Aalto et al. 1991a; Garay et al. 1993) 2 . 
We measured total 12 CO(2-l) and 13 CO(2-l) flux of 4.1 xlO 3 Jy km s' 1 and 1.4xl0 2 Jy 
km s _1 from our natural-weighting data corrected for the primary beam attenuation and 
convolved to the SEST resolution, which we assumed to be 23" in FWHM at 230 GHz. The 
high fraction of flux recovery is consistent with the compact size of the CO emission, whose 
half-power diameter is estimated to be about 10" by Aalto et al. (1991a). The 20% larger 
flux in 12 CO at the SMA than at the SEST may be due to a larger-than- usual error in the 
SMA calibration for this source at very low elevations, and may be also due to pointing and 
calibration uncertainties in the single-dish observations 3 . The single-dish observations must 
also have slightly underestimated the CO flux by over-subtracting the emission as a linear 
baseline. This is because the CO line is wider than previously believed as we see below. How- 
ever, this leads to only 1-2% underestimation of the single-dish CO flux. In any case, our 
relative calibration between the CO lines and continuum should be more accurate than their 
absolute calibration, because they were simultaneously observed with the same receivers and 
went through the same signal path and reduction procedure. Regarding missing flux in our 
higher resolution datasets, our highest-resolution 12 CO data cube contains 54% of the total 
flux in our natural-weighting data. 



3. Overall Gas Distribution and Kinematics 

3.1. spatial distribution of molecular gas and dust 

The 12 CO (2-1) emission as shown in Fig. 2 is highly concentrated toward the galactic 
center. It has a half-maximum radius of about 5" (~1 kpc) in the galactic plane (Fig. 3). 
The compact CO morphology is consistent with the SEST mapping result of Aalto et al. 
(1991a). The CO central peak is elongated in the north-south direction and connects the 
double nucleus of 5" separation. There are two peaks in this structure in the 13 CO(2-l) map 
also shown in Fig. 2. Each 13 CO peak is spatially associated with one of the two nuclei. We 
discuss these double peaks in more detail in the following sections using higher resolution 
images. The gas distribution extends at least to the radius of 3 kpc in our data, and it may 



2 Only the 12 CO(2-l) flux of Aalto et al. (1991a) is cited here because their 13 CO(2-l) flux is more than 
twice larger than the other two observations. 

3 The size of the SEST beam given in the papers ranges from 22" at 1.3 mm (=230.6 GHz) to 25" at 230 
GHz. The conversion factor from temperature to flux density scales with the square of the beam size for a 
compact source. Thus the single-dish flux values given above increase by 18% if we adopt the largest beam 
size. 
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be extended beyond our field of view. In the outer area, there are arcs or hints of spiral 
arms in the molecular gas. Some of them correspond to spiral arms or dust lanes in the 
optical images in the literature. For example, the horizontal feature near the northeast edge 
of the 12 CO map, around the offset of AR.A.=10" and ADec.=+15", corresponds to a spiral 
arm emanating from the galactic center (see the Ha image of Moorwood & Oliva 1994, Fig. 
2(b)). 

The 1.3 mm continuum and C 18 0(2-l) emission are also detected in the galactic center 
(see Fig. 2). The spatial distribution of the continuum is quite similar to that of 13 CO(2-l) 
emission. There are peaks at or near the two nuclei, surrounded by extended emission. A 
weaker peak is also seen ~6" east of the nucleus N, as is seen in 13 CO (2-1) and 12 CO (2-1). 
The shape of the C 18 0(2-l) emission looks somewhat different from those of other emission 
mentioned above. However, this is likely due at least partly to the low signal-to-noise ratio 
in the data. 

The total 12 CO(2-l) fluxes within the galactocentric radii (r g ) of 10" and 20" are 3.4x 10 3 
Jy km s _1 and 5.2 x 10 3 Jy km s -1 , respectively, as shown in Fig. 3. The total fluxes of 
13 CO (2-1) and C 18 0(2-1) emission within the galactocentric radius of 10" are 1.5 x 10 2 
Jy km s _1 and 31 Jy km s _1 , respectively. The total flux density of the 1.3 mm (222.7 
GHz) continuum in the same area is 79 mJy Each line flux is integrated over the same 600 
km s _1 range from 2530 km s _1 to 3130 km s _1 , and each flux measurement for both line 
and continuum was made in the same way as explained in the caption of Fig. 3. Each flux 
has a 10% scaling uncertainty as noted in §2. 



3.2. gas mass estimate 

The mass of molecular gas (M mo i) in the region is estimated in various ways from the 
line emission. The molecular gas mass of 2 x 10 10 M & for r g < 10" = 1.7 kpc is obtained from 
the conversion factor Nn 2 / Ico(i-o) — 3.0 x 10 20 cm -2 (K km s" 1 ) -1 , which is based on the 
virial analysis of molecular clouds in the disk of our Galaxy (Scoville et al. 1987; Solomon 
et al. 1987). We assume that the CO(1-0) and CO(2-l) lines have the same brightness 
temperature and emitting area. This mass estimate and the ones in the rest of this section 
include the contribution of helium. The M mol in the same region is estimated to be lower, 
9 x 10 9 M Q and 2 x 10 9 M , from the conversion factors derived from 7-ray observations 
to count H-nuclei; the two estimates are from the conversion factor of 1.56 x 10 20 cm~ 2 
(K km s -1 ) -1 assumed to be constant over the Galaxy (Hunter et al. 1997) and 0.4 x 10 20 
cm -2 (K km s -1 )^ 1 at r g = 2 kpc in our Galaxy (Strong et al. 2004), respectively. Another 
estimate is from the C 18 0(2-l) emission and the prescription of Stutzki & Glisten (1990) 
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and Wild et al. (1992). They suggested that this optically thin line should have nearly 
a constant emissivity per molecule in the range of temperature and density relevant for 
starburst galaxies such as M82. NGC 3256, also being a starburst, is probably in the class 
of galactic nuclei to which the method is applicable. The observed C 18 0(2-l) flux gives the 
molecular gas mass of the 1.4 x 10 9 M Q for r g < 1.7 kpc for the assumed C 18 abundance 
of [C 18 0]/[H2] = 10~ 6 ' 8 (Frerking et al. 1982). The mass of atomic gas within this radius 
is estimated to be M atom =(0.4-1.9) x 10 9 M by scaling the mass estimate by English et al. 
(2003), who obtained the mass of HI absorbing gas in the central 23". 

The 1.3 mm continuum is also used for the mass estimate in the following way. First, 
about a half of the 1.3 mm flux density, or 45 mJy, is attributed to dust emission by sub- 
tracting the contributions from synchrotron and free-free emission from the observed flux 
density. The synchrotron emission at 1.3 mm is estimated to be 15 mJy by extrapolating 
a power law (of a spectral index of —0.76) obtained by fitting the radio flux measurements 
between 160 MHz and 5 GHz in the NASA Extragalactic Database 4 . The free-free flux den- 
sity is estimated to be 19 mJy from the H II region model of Roy et al. (2005), who derived 
properties of the H II regions in the galactic center by observing radio recombination lines in 
NGC 3256. After subtracting these contributions, the remaining emission of 45 mJy must be 
the thermal emission from dust in the interstellar medium. We used the dust emissivity of 
Draine & Lee (1984) and the dust temperature of 44 K (Smith & Harvey 1996) to obtain the 
gas (molecular plus atomic) mass of 2.6 x 10 9 M Q for r g < 1.7 kpc. The mass is reduced by 
half if we adopt the dust opacity coefficient of Hildebrand (1983). An assumption made in 
these estimates is that the gas to dust ratio and dust properties are the same in our Galaxy 
and NGC 3256. 

The large scatter among the mass estimates is not surprising for a starburst nucleus. 
Many assumptions that can not be easily verified are inevitably made, including those on 
metallicity, molecular abundance, dust properties, and the properties of the molecular gas. 
For example, the high pressure in the interstellar medium as implied by the superwind from 
the starburst probably helps to bind molecular clouds. This increases the 12 CO emissivity 
of the clouds and leads to an overestimation of gas mass if the conversion factor derived 
for selfgravitationally bound clouds is applied (Bryant & Scoville 1996; Oka et al. 1998). 
Independent of these assumptions, dynamical mass sets a stringent constraint since the gas 
mass can not exceed it. As we see in the next section (§3.3), the dynamical mass for the 
central region is estimated to be ~1.2 x 10 10 M Q for r g < 1.7 kpc. 



4 The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, Cali- 
fornia Institute of Technology, under contract with the National Aeronautics and Space Administration. 
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We adopt the total mass of neutral gas in the central region of M gas (r g < 1.7 kpc) ~ 
2.9 x 10 9 M Q which is the logarithmic average of the three mass estimates from 12 CO(2-l), 
C 18 0(2-l), and dust emission. For each of the three methods, we geometrically averaged 
the estimates given above, except the one based on 12 CO that exceeded the dynamical mass. 
We assign our estimate a factor of 2 error in each direction, on the basis of the range of the 
three estimates. 



3.3. kinematics of the main gas disk 

The velocity field of the 12 CO(2-l) emission in Fig. 2 shows a spider pattern that is 
typical of gas disks in spiral galaxies. The kinematical major axis is at a position angle of 
~90° in the vicinity of the double nucleus, and appears to be ~70° in the outer regions. 
The outer position angle is consistent with those from previous spectroscopic observations 
in CO(2-l) (Aalto et al. 1991a) and Br7 (Moorwood & Oliva 1994). It also agrees with the 
isophotal fitting in the -fT-band by Rothberg & Joseph (2004), in which the major axis has 
a nearly constant position angle of ~60° between the sky-plane radius of 2 kpc and 6 kpc. 
The overall sense of rotation, i.e., east side being redshifted and west blueshifted, is the same 
as that of H I tidal tails (English et al. 2003). The apparent ~20° shift in the kinematical 
major axis could be due to streaming motion of the gas or warp of the main gas disk. The 
velocity field seen at this spatial scale (i.e., kpc scale) is relatively undistorted for a merging 
galaxy. 

The cold gas component appears to have kinematically merged (or relaxed) and formed 
a disk in the kpc scale. We therefore fitted the 12 CO velocity map in Fig. 2 using the 
AIPS task GAL to derive the kinematical parameters of the gas disk; the accompanying 
intensity map was used for data weighting. The intensity weighting makes our parameters 
more appropriate for the inner part (r < 10") of the CO disk, reducing the effect of possible 
warp or non-circular motion in the outer disk. The resulting dynamical center is between 
the two nuclei N and S. In order to assess the robustness of the fitting results by altering 
the effect of beam smearing, we did the same fitting using the moment maps made with the 
uniform weighting, which gave a better spatial resolution of 3'/3 x l'/3 while retaining 87% 
of the flux in the natural weighting (4'/8 x l'/9) data. The dynamical center, inclination, 
position angle of the major axis, and the systemic velocity agreed within 0'/8, 8°, 2°, and 
9 km s _1 , respectively. The average parameters are ckm — 10 h 27 m 51!25, 5m = — 43°54'15 / .'7 
(J2000), im = 45°, P.Am = 94°, and V^ ySj M = 2791 km s -1 , respectively, where the subscript 
M is for the main (or merged) disk. One half of the abovementioned difference between the 
two fitting results gives each average value an estimate of uncertainty, though it is probably 
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underestimated because any systematic effect such as non-circular or non-coplanar motion 
of gas could bias our fitting results. 

The parameters we estimated are in reasonable agreement with those from previous 
observations. Specifically, the inclination and position angle are consistent with those by 
Feast et al. (1978) and others. The systemic velocity is consistent with those by Casoli et 
al. (1991), Aalto et al. (1991a), and many others. We use in this paper V sys = 2775 km s _1 
taken from the literature. No more than 10% error is introduced in the velocity related 
parameters derived in the following — such as dynamical mass, momentum, and energy - 
by the slight difference (16 km s _1 ) from our fitted value. Lastly, but most notably, the 
dynamical center between the two nuclei is consistent with that measured from Ha velocity 
field by Lipari et al. (2000). The parameters we derived from the millimeter- wave CO line 
are not affected nor biased by dust obscuration, which is known to be large in the southern 
part of the main disk. (For example, Kotilainen et al. (1996) estimated Ay ~ 10.7 mag 
for line emission toward the S nucleus.) As cautioned, the decoupled dynamical systems we 
suggest to exist around the double nucleus (see below) may have affected the parameters to 
some extent. We expect this effect to be small, because local distortion of the velocity field 
is not obvious around the double nucleus in the low resolution maps. 

The dynamical mass within the radius of 10" =1.7 kpc of the main gas disk is estimated 
to be M dyn (M, r < 1.7 kpc) = (1.2 ± 0.4) x 10 10 M from the rotation velocity V(r = 
10") suum — 125 ± 20 km s" 1 and the inclination (45°± 4°) obtained from the kinematical 
fit. (See Fig. 6 (c) for the rotation velocity.) The fraction of gas in the dynamical mass is 
in the range of 0.1-0.5 for our estimate of gas mass and the central value of the dynamical 
mass. 



4. Double Nucleus 

4.1. gas distribution around the double nucleus 

The high-resolution map in Fig. 4 reveals compact gas peaks near or around the nuclei 
N and S. The southern peak is more compact than the northern one, which is more extended 
than the synthesized beam. The two gas peaks are connected with a bridge of CO emission. 
There are also arcs or partial spiral arms as well as secondary peaks that were already hinted 
in the lower resolution maps. The CO distribution traces some of the conspicuous dust lanes 
in the optical, as shown in Figure 5 that compares the CO contours with the /-band HST 
image of Zepf et al. (1999). Most notably, the following features correspond to dark lanes 
in the HST image: the arc-like dust lane about 4"-8" east of the northern nucleus, the gas 
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peak about 5" east of the southern nucleus, and the gas peak around the southern nucleus 
itself as well as the westward tail emanating from the peak. The 'tail' feature of molecular 
gas actually trails the S nucleus, judging from the clockwise rotation of the galaxy inferred 
from the spiral patterns in the HST image. 

The 12 CO (2-1) integrated intensity averaged within a 2"-diameter aperture centered on 
each nucleus is 44 Jy km s" 1 arcsec" 2 for both of the nuclei. This corresponds to a surface 
density of neutral gas 1.3 x 10 3 M pc~ 2 on the sky plane if we use the same scale as our 
mass estimate in the previous section and ignore the correction for missing flux. It also 
corresponds to a visual extinction of Ay ~ 60 mag for a uniform slab in front of the stars. 
The actual extinction toward the nuclei is obviously much smaller than this because at least 
near-IR images show both nuclei. Kotilainen et al. (1996) estimated A v ~ 2.4 mag and 
~ 5.3 mag for the continuum emission from N and S nucleus, respectively. The reasons for 
these smaller extinctions are most likely that stars and gas are mixed, that the gas and dust 
has a smaller scale height than the stars, and that the molecular gas consists of clumps. 

It is notable, however, that the S nucleus suffers more extinction than the N nucleus 
despite the same gas surface densities inferred in the directions of the two nuclei. The 
absorbing column density suggested from the actual extinction is a factor of 2.2 larger for 
the S nucleus than for the N nucleus. This discrepancy suggests a geometrical effect. The S 
nucleus is probably more deeply embedded in gas and dust or located behind the main gas 
disk. 

The majority of 12 CO(2-l) emission in the galactic center comes from the main gas 
disk, despite the peaks around the two nuclei in the integrated intensity map. The areas 
within 1" from the two nuclei have only 8 % of the 12 CO(2-l) flux observed in the r g < 10" 
region, which is the area within 10" = 1.7 kpc from the midpoint of the double nucleus in 
the plane of the main gas disk. This fraction doubles, but is still small, if we apply a uniform 
correction for the flux resolved out in the highest resolution map. 

4.2. gas kinematics and dynamical masses 

4-2.1. disturbances in the main disk 

The high-resolution mean velocity map in Fig. 4 shows the same overall rotation of the 
main gas disk as the lower resolution data. In addition it reveals small scale disturbances 
in the disk. The most notable of these outside the double nucleus is the noncircular motion 
associated with the spiral feature about 4"-8" east of the N nucleus. The merged gas disk 
is certainly disturbed at the few 100 pc scale. However, we do not see in the main disk a 



- 12 - 



double-peaked line profile that would have suggested two distinct components overlapping in 
our line of sight. It appears that most of the gas from the merger progenitors has coalesced 
in the main gas disk, except perhaps in close proximity to the two nuclei. 

4-2.2. possible mini-disks around the double nucleus 

On, and very close to, each nucleus, the velocity map indicates a steep velocity gradient 
in the east-west direction across the nucleus. The steep velocity gradients are confirmed in 
the position-velocity (PV) diagrams shown in Fig. 6. The full velocity width within 1" of 
each nucleus is about 250 km s" 1 and 310 km s _1 for the N and S nucleus, respectively. Our 
data do not show a clear sign of velocity offsets, i.e., relative motion, between the main gas 
disk and each nucleus. 

We model the velocity gradient as rotation of gas around each nucleus. The Keplerian 
dynamical mass within 1" (= 170 pc) of each nucleus is Md yn (N, r < 170 pc) = 6 x 10 8 sin -2 i^ 
M & and Md yn (S, r < 170 pc) = 9 x 10 8 sin~ 2 «s M@ in this model, where i^ and is are the 
inclinations of the spin axes with respect to our line of sight. The uncertainty in these 
dynamical masses is probably as large as a factor of 2 reflecting the uncertainty in the radii 
of the peak velocities. For i N ~ i s ~ 45°, the dynamical masses are approximately 1 x 10 9 
M and 2 x 10 9 M for N and S, respectively. 

The dynamical mass for the N nucleus is consistent with (i.e., reasonably larger than) 
the optical estimate, ~10 8 M Q within 40 pc, based on the HST STIS spectra across the 
nucleus (Neff et al. 2003; Lipari et al. 2004). The mass of neutral gas within 1" radius in the 
suggested disk is about 1 x 10 8 cosi N M for the N disk (and 1 x 10 8 cosi s M Q for the S) 
if we scale the gas mass estimated in §3.2 with 12 CO (2-1) flux. The gas-to-dynamical mass 
ratio therefore does not prohibit the mini-disk model. The total mass of the two nuclei is at 
least 15 % of the dynamical mass within a 10" (1.7 kpc) radius from the dynamical center 
of the merger. 

In the mini-disk model, the suggested rotation axis of each nucleus is roughly parallel 
to that of the main gas disk, and the sense of velocity gradient is about the same in the 
three disks with eastward being redshifted and westward blueshifted. This is consistent 
with the prograde-prograde merger that has been suggested for the galaxy. In the merger 
configuration, the spin axes of the progenitors and their orbital rotation axis are roughly 
parallel with each other and all rotation is in the same direction. Thus the dense gas disk 
to form around the core of each progenitor and the large gas disk to form from merged disk 
gas will have spin vectors in roughly the same direction and orientation, as in NGC 3256, if 
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the transfer of angular momentum is negligible between the spin and orbital motion of the 
mini disks. 

Obviously one has to be careful about the interpretation of the velocity structure in 
the main gas disk because of the unfortunate configuration that the two nuclei are aligned 
along the minor axis of the main disk. This configuration makes it difficult to separate the 
rotation of the main disk and that around each nucleus for a prograde-prograde system. We 
therefore made a position-velocity diagram across the midpoint M of the N and S nuclei. 
The PV diagram shown in Fig. 6 (c) still shows the steep velocity gradient that we saw on 
the nuclei N and S. There is, however, also a component with a shallower velocity gradient 
that almost linearly rises to the velocity offset of ~120 km s _1 at the offset of 10". We infer 
that this component with the shallow velocity gradient represents the main gas disk and 
that the central component with a steep velocity gradient is the contamination from the N 
and S nuclei. The contamination is quite possible because the midpoint M is only 2'.'5 from 
each nucleus while our spatial resolution in the north-south direction is 2". 

This interpretation is supported by the PV diagram through the two nuclei and their 
midpoint (Fig. 6 (d)). It suggests wider line widths at the two nuclei than at their midpoint. 
Additional supporting evidence for the shallow rise of the rotation curve of the main disk is 
that the isovelocity contours in the mean velocity map (Figs. 2 and 4) are roughly parallel 
to the minor axis within ~5" from the dynamical center of the main disk. This indicates 
a near rigid-body rotation of the main disk in the region. It is consistent with the shallow 
velocity gradient in Fig. 6 (c) but not with a steeply rising rotation curve that approaches 
near its maximum velocity at a radius of ~1". 

Another cautionary remark is due for the outflow of ionized gas that Lipari et al. (2000) 
suggested to emanate from the N nucleus. We regard it as unlikely to be the cause of the steep 
velocity gradient of molecular gas across the nucleus. This is because the reported outflow is 
several kpc large and has an outflow axis in the position angle of 150° - 160°. The flow axis 
is not aligned with the velocity gradient of the molecular gas but roughly perpendicular to 
it. Thus, although a starburst-driven wind can entrain molecular gas, the one reported from 
the optical observations would not cause the observed CO velocity gradient. Our estimate 
of the dynamical mass of the nuclei is therefore not affected by the outflow. No report of an 
outflow has been made for the S nucleus so far, although, if it exists, it would not be easily 
detected through the large extinction. 

A qualitative argument supports the suggested gas disks around the two nuclei, though 
observational confirmation of them requires data at a higher resolution. If near-IR and radio 
peaks N and S are indeed the remnant nuclei of the merger progenitors that currently orbit 
in (or in and out of) the main gas disk, then collisional gas can follow each stellar nucleus 
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to make a peak around it only if the gas is gravitationally bound to the nucleus. A massive 
remnant nucleus sinks toward the dynamical center because of dynamical friction. The gas 
around the nucleus tend to be left behind if not bound to the nucleus, because gas clouds 
are less massive than the nucleus and also because of hydrodynamical friction between the 
gas around the nucleus and that in the main gas disk. For cold molecular gas bound to each 
nucleus, the stable configuration is a rotationally supported disk around the nucleus. In this 
picture, some of the gas around each nucleus can be captured from the main gas disk. 

4-2.3. comparison with other galaxies 

Small gas disks around merger nuclei separated by less than 1 kpc have been observed 
in Arp 220 (Sakamoto et al. 1999; Mundell et al. 2001). Such a mini disk was also suggested 
in the center of M83 and attributed to a minor merger (Sakamoto et al. 2004). In Arp 220, 
the dynamical mass of each disk is at least 1-2 x 10 9 M Q within a radius of 100 pc and the 
two nuclei are separated by 0.3 kpc on the sky. 

Another case of some similarity with NGC 3256 is the luminous infrared galaxy NGC 
6240, which is a merger with two nuclei separated by 750 pc on the sky. Tecza et al (2000) 
found a steep velocity gradient across each nucleus in the stellar velocity field measured with 
a CO absorption band, and inferred a pair of rotating stellar systems at the two nuclei with 
a dynamical mass for each of 2-8 xlO 9 M within a radius of 200 pc. 

The pairs of compact mass concentrations in the two major mergers Arp 220 and NGC 
6240 are comparable in size and mass with that in NGC 3256. NGC 3256 may be a new 
member of the (yet) small group of merging galaxies that show massive and kinematically- 
decoupled rotating disks in their central kpc. 

4.2.4. alternative models 

Lastly, there could be a totally different interpretation of our high-resolution observa- 
tions. Namely, one could assume that the putative nuclei N and S are not really remnant 
nuclei of the progenitors but just bright star forming regions and that they are there because 
there are compact concentrations of gas for star formation (c.f., Eckart & Downes 2001). 
This alternative would need a model to explain why there are two gas peaks in this galaxy 
and why there appears to be large velocity gradients across the gas clumps. We have not 
attempted such modeling. If both N and S have a low-luminosity AGN as suggested by Neff 
et al. (2003), then they are most likely galactic nuclei. 
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We also did not explore the opposite model by Lipari et al. (2000) that NGC 3256 
has a third merger nucleus at 6" east of the N nucleus. The location is close to the CO 
spiral feature that we noted with its non-circular motion. However, the non-circular motion 
appears to be associated with the spiral rather than localized on the suggested nucleus. The 
line width at the position, at —6" in Fig. 6 (a), is less than half of the ones at the N and S 
nuclei. It is not significantly larger than the ~100 km s _1 line widths seen at other locations 
in the turbulent main gas disk (see Figs. 6 (b) and (c)). The third nucleus, if exists, appears 
to play a minor role in gas dynamics. 



5. High- velocity gas 

5.1. observational signatures 

We have discovered high velocity molecular gas in the galactic center. This is most 
clearly seen in the position-velocity diagram in Figure 7 as the high velocity component at 
the offset of 0". The high velocity gas is seen in both blueshifted and redshifted velocities 
and is more prominent in the latter. The 12 CO PV diagram is along the major axis of the 
galaxy and is made by integrating the 6" wide area along the line of nodes. The double 
nucleus is within this 6"-wide 'slit' and is also at the offset of 0". The data cube used for 
this PV diagram has 4'/7 x l'/8 (FWHM) and 10 km s -1 resolutions and has 92 % of the 
12 CO flux in the natural- weighting dataset. In this dataset, the redshifted gas is detected up 
to the velocity of 3120 km s _1 , which is about 350 km s -1 from the systemic velocity of the 
galaxy. The blueshifted gas at the galactic center is detected down to the velocity of ~2580 
km s _1 , or about 200 km s -1 from the systemic velocity, though an arm-like feature about 
15" west of the double nucleus has velocities down to about 2540 km s^ 1 . We searched in 
the range of 2400 - 3320 km s" 1 for emission at even more extreme velocities, but did not 
detect any at the abovementioned resolutions. The redshifted component, however, appears 
to have a faint tail up to the velocity of ~3200 km s _1 in a lower-resolution dataset as we 
see below. We did not detect the high velocity gas in 13 CO and C 18 0. 

The high velocity gas is located between the two nuclei. This is seen in Figure 8 that 
shows 12 CO(2-l) maps made by averaging the velocities of the redshifted and blueshifted high 
velocity gas. In both maps, the high velocity gas is not peaked on either of the two nuclei but 
between them. The peak signal-to-noise ratio of the high velocity gas in the maps is 12.8 and 
4.4 for the redshifted and blueshifted gas, respectively. Any error in continuum subtraction 
can hardly affect this emission, because the maximum flux density of the continuum is less 
than a third of the peak flux density of the fainter, blueshifted high velocity gas. The peak 
position of the redshifted high velocity gas is a = 10 h 27 m 5r!22 and 5 = -43°54'16'/3 (J2000) 
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according to a Gaussian fit. The emission has a deconvolved size of 4'/6 x 2'/5 (P.A.=161°). 
The emission is too weak, and possibly too extended, to be detected in our uniform weighing 
maps used in the previous section. 

The high velocity gas appears almost at the same position in the redshifted and blueshifted 
velocities. Although the blueshifted emission is too weak to precisely determine its centroid 
position, the peak positions in the two maps in Figure 8 agree within l'/3. In addition, we 
did not find a systematic velocity gradient within the redshifted high velocity gas. To check 
this, we made a mean velocity map of the high velocity gas to see velocity gradient in any 
particular direction, and also inspected channel maps to see if there was a systematic shift 
of the emission centroid. We did not detect a sign of velocity gradient in either of them. 
This can constrain the parameters in models for the high velocity gas (e.g., rotation curve 
for a rotation model and outflow geometry for an outflow model). 

The integrated 12 CO(2-l) flux of the high velocity gas is 3 % of the total 12 CO(2-l) 
flux in the 5"-diameter area at the galactic center. This is measured from the spectrum in 
Figure 9. The spectrum is extracted at the midpoint of the double nucleus from our natural 
weighted data convolved to the 5" resolution. The emission line spans from 2530 km s _1 to 
3200 km s _1 when we use the first null on each side to determine the line edge. The full 
width at zero intensity (FWZI) of the emission is thus 670 km s _1 . The CO flux in the red 
component above 2990 km s _1 is 23 Jy km s _1 in the central 5", that in the blue component 
below 2610 km s _1 is 5 Jy km s _1 , and the total flux in the line is 936 Jy km s _1 . 

5.2. molecular outflow model 

5.2.1. supporting evidence for outflow 

We suggest that the high velocity g bipolar outflow of molecular gas from the 

main gas disk. This is mainly because the galaxy is known to have a superwind extending 
several kpc out of the central starburst (Scarrott et al. 1996; Moran et al. 1999; Heckman et 
al. 2000; Lipari et al. 2000) and the flow of hot gas can entrain cold molecular gas. It is also 
because alternative explanations of the high velocity gas seem less likely as we see below. 
The location of the high velocity gas makes the main gas disk the most likely place from 
which the molecular wind emanates. 

Among the existing observations of the superwind, Heckman et al. (2000) detected to- 
ward the central few arcsecond of NGC 3256 a doublet of Na I D absorption line blueshifted 
by about 300 km s~ x with respect to the systemic velocity of the galaxy. In our CO observa- 
tions, the maximum velocity offsets are +420 and —240 km s -1 with respect to the systemic 
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velocity, according to the FWZI line width. The velocity offsets of the Na I and CO lines are 
in good agreement in magnitude. Both lines trace neutral gas; the former traces atomic gas 
and the latter molecular gas. Thus the agreement supports the outflow model of the high 
velocity molecular gas. 

5. 2. 2. alternatives 

Before going further into the outflow model, we briefly check its alternatives. First, 
rotation can cause high velocity. It would have to be the one around the dynamical center of 
the merger, where the high velocity gas is. The dynamical mass needed for the high velocity 
would be at least 3 x 10 9 M & within a radius of 110 pc (=0'.'65) if we use the CO line width 
for velocity and the nominal offset between the blueshifted and redshifted CO peaks for size. 
The dynamical mass is larger than those of the two nuclei (§4.2.2) despite the lack of a 
X-band peak between them. In addition, the large mass does not fit the rotation curve that 
we inferred in §4.2.2 to be nearly linearly rising to the radius of 10". The dynamical mass 
can be smaller and consistent with the observations if the extent of the high velocity gas is 
smaller, but a smaller gas disk means a higher mass density and a higher CO intensity in 
it. For example, a disk with a 10 pc radius would require the dynamical mass, mass density, 
and brightness temperature of 3 x 10 8 M & , 5 x 10 4 M pc~ 3 , and 100 K, respectively. An 
AGN may realize these conditions, but none has been found between the two nuclei. Thus 
we regard the rotation model less likely than the outflow model. 

Other alternatives include locally large velocity dispersion, gas infall from out of the 
main disk, and an outflow driven by an AGN jet. The velocity dispersion does not seem 
to have a local source of turbulence, unless there is a special mechanism related to the 
merger gas dynamics. The gas infall from both (i.e., approaching and receding) sides toward 
the dynamical center seems too much a coincidence. No AGN jet has been seen in the 
region. Note that all the alternative models except the jet-driven outflow can not explain 
the blueshifted Na I absorption, though the CO and Na I features could be of different 
origins. 

5.2.3. outflow parameters 

In the outflow model, we can estimate the following parameters from our observations. 

Velocity: - If the outflow is well collimated along the rotation axis of the main gas 
disk then the maximum flow velocity corrected for inclination would be 600 km s _1 for the 
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receding flow and 350 km s 1 for the approaching flow. However the outflow may have a wide 
opening angle. The inclination correction in that case would be smaller or even negligible. 

Mass: — The mass of molecular gas involved in the high-velocity flow is computed to 
be 1.4 x 10 7 M Q and 0.3 x 10 7 M for the receding and approaching gas, respectively, from 
the 12 CO(2-l) flux in the velocity ranges where the high velocity gas is distinguishable from 
the disk gas. These estimates are based on the 12 CO(2-l)-to-gas mass scaling factor that 
we adopted for the galactic center (§3.2). Their uncertainties are larger than that for the 
disk gas because we do not know whether and how the 12 CO (2-1) emissivity changes in the 
outflow gas. The contribution from atomic gas, which was a factor of 0.3 in the gas mass 
estimated in §3.2, is excluded from the above numbers because atomic-to-molecular mass 
ratio in the outflow is likely different from that in the disk gas. The numbers given here 
suggest that molecular gas of the order q/10 7 M Q is being ejected from the central kiloparsec 
at velocities > 200 km s _1 in the molecular outflow. There is likely a larger amount of 
molecular gas in the outflow at smaller velocities. 

Momentum and Kinetic Energy: — The momentum and kinetic energy of the molecular 
outflow are estimated using the same mass scaling and integration on the spectrum in Fig. 
9. The red and blue high velocity gas have momentum of 8 x 10 42 kg m s -1 and 1 x 10 42 kg 
m s _1 , respectively, in the direction of our line of sight. The kinetic energy of the redshifted 
high-velocity gas is 1 x 10 48 J, without correction for the inclination of the flow, while the 
blueshifted gas has one tenth of that energy. The unknown contribution of outflow gas slower 
than 200 km s -1 becomes smaller for the momentum and kinetic energy because of the (Av) 1 
and (Av) 2 weights used for these moments. 

Timescale : — We estimate the characteristic timescale for the redshifted outflow with 
more than 200 km s" 1 velocity offset to be ~2 Myr from the deconvolved size of the high 
velocity emission and the outflow velocity. This is not the age of the outflow, which should 
be comparable to the age of the starburst and is probably an order of magnitude larger. The 
timescale is rather the crossing time, for the high- velocity gas, of the < 500 pc region where 
the high velocity emission is detectable. For symmetry, we assume the same timescale for 
the blueshifted high velocity gas. 

Fluxes and Mechanical Luminosity: — We can estimate the mass flux, momentum flux, 
and mechanical luminosity of the high velocity molecular outflow using the timescale and 
the parameters estimated above. The mass flux of molecular gas in the high velocity outflow 
is ~9 M Q yr -1 , the momentum flux injected to the outflow is ~1 x 10 29 N, and the kinetic 
luminosity of the outflow gas is ~2 x 10 34 W. These values are for the high velocity molecular 
component. Contributions from atomic and ionized gas and those from slower velocity flow 
increase these numbers. For comparison, the average rates of 17±f M Q yr" 1 , 10 29 - 1±0 - 5 N, 
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and io 346±0 - 6 W have been reported for a sample of infrared luminous galaxies (io 1L36±0 - 4 
I/© ) from Na I observations (Rupke et al. 2005). 

The uncertainties in the parameters derived here are worth summarizing. The param- 
eters involving gas mass inherit the uncertainties of the mass estimate, with an additional 
uncertainty from the possible variation of gas properties between the disk and the outflow. 
Those involving velocity have uncertainty due to the unknown outflow geometry. We did not 
apply an inclination correction, in effect assuming a wide opening angle of the flow. Still, 
velocity-related parameters are estimated from velocity components along our line of sight. 
If we apply the correction for a well-collimated outflow 45° inclined to our line of sight, then 
momentum and kinetic energy increase by a factor of 1.4 and 2, respectively. The momentum 
flux and the kinetic luminosity also increase by a factor of 1.4 and 2, respectively, since the 
crossing time does not change with the correction. As noted, the parameters do not include 
the contribution of the outflow gas whose low line-of-sight velocity does not distinguish the 
gas from the disk gas. 



5.2.4- energy and gas consumption budgets 

The energy of the outflow can be supplied by a portion of the starburst. If we equally 
divide the far-IR luminosity of the galaxy among each of the two nuclei and the main disk, 
then the mechanical luminosity released from the starburst in the main disk is ~3 x 10 35 
W assuming a 10 Myr-old continuous starburst that has the Salpeter initial mass function 
(IMF) in the mass range of 1-100 M Q (Leitherer et al. 1999). Thus the molecular outflow is 
energetically possible with a wide margin. It can be caused by passing < 10 percent of the 
mechanical energy from the disk starburst to molecular gas; the rest can be used to drive 
the superwind of hot gas. 

The star formation rate needed to generate the starburst luminosity is ~10 M yr _1 
under the same assumptions on the starburst. The rates are therefore of the same order 
of magnitude for the consumption of molecular gas by star formation and for the dispersal 
of molecular gas by the outflow. This comparison, however, should be viewed with caution 
because the star formation rate estimated from luminosity, as well as the gas outflow rate, has 
a large uncertainty. The former rate can change by a factor of 10 with different assumptions 
about the IMF, because high mass stars generate most luminosity while low mass stars 
usually determine most of the gas consumption. 

Perhaps a more robust number to describe the significance of the molecular outflow is 
the depletion time of molecular gas solely by the outflow. It is estimated to be 70 (= 2/0.03) 
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Myr from the fraction of high velocity gas in the central 5", 0.03, and the time scale of 
the outflow, 2 Myr, under the simplistic assumption that the high velocity gas does not 
return 5 . The timescale is comparable to the typical timescale of a starburst, again implying 
the significant role of the molecular outflow in the starburst. This timescale of gas depletion 
or dispersal by the outflow is estimated without explicitly using the CO-to-H2 conversion 
factor. Hence the uncertainty in the depletion time due to the conversion factor is only 
indirect, through any variation of the scaling factor between the outflow and disk gas. We 
caution that the timescale is derived to measure the significance of the molecular outflow 
rather than to predict the evolution of the system, for which one needs to consider such 
factors as the time evolution of star formation, that of the outflow, and further gas accretion 
within the galaxy. 

To summarize, our observations of the molecular outflow agree with the optical absorption- 
line surveys toward starburst galaxies (Heckman et al. 2000; Rupke et al. 2005) in that a 
starburst-driven superwind can be as important as the star formation itself in the gas con- 
sumption budget of a luminous infrared galaxy. Rupke et al. (2005) obtained a median 
value of 0.33 for the mass entrainment efficiency (i.e., the rate of mass loss due to super- 
wind normalized by the star formation rate) in luminous infrared galaxies, to which NGC 
3256 belongs. The efficiency for the molecular outflow in NGC 3256 is ~1 if a third of the 
starburst in the galaxy contributes to the observed outflow, or ~0.3 if the entire starburst 
contributes. Thus the entrainment efficiency for molecular gas in NGC 3256 agrees, in the 
order of magnitude, with the median efficiency estimated for atomic gas. 



5.2.5. geometry 

We finally comment on the spatial configuration of the starburst in the galactic center. 
The spatial distribution of massive star formation is peaked toward the N nucleus as seen in 
[Fe II] and Br7 (Moorwood & Oliva 1994), Paa (Alonso-Herrero et al. 2002), Ha and diffuse 
X-ray in the 0.3-10 keV range (Lira et al. 2002). The S nucleus also shows a peak in Br7, 
though its flux corrected for extinction is estimated to be 4 times smaller than that of the 
N nucleus (Kotilainen et al. 1996). In contrast, the region between the two nuclei does not 
have a peak in these tracers of massive star formation, even though there are H II regions 



5 We can not easily assess whether the molecular gas escapes the galaxy or not, without information 
about the current height of molecular gas above the main disk, the shape of the gravitational potential, 
and structure of the superwind. If the gas returns, then the timescale is the one during which all the gas 
experiences the out-of-plane trip at least once. The molecular gas may well be dissociated and ionized in 
the hot wind during the trip. 
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and super star clusters spreading between and around the two nuclei (Alonso-Herrero et 
al. 2002). It is therefore somewhat puzzling to find a molecular outflow between the two 
nuclei and not from each nucleus, in particular from the nucleus N. Unfortunately the optical 
observations that indicated the presence of a superwind do not tell us the exact geometry of 
the wind and whether there are multiple superwinds in the region. 

Clues to solve the puzzle may be that the base of the superwind likely has a size of 
1 kpc or larger and that the two nuclei as well as the star forming clusters in the region 
have been moving with respect to the dynamical center of the system. Heckman et al. 
(1990) measured the pressure profile in the center of NGC 3256. They found that the 
central 1-2 kpc diameter region has very high pressure (~10 -9 ' 5 Pa) with a shallow decline 
outwards, and that the pressure gradient steepens outside the region. As they interpreted 
the profile, this is expected for a superwind with a base diameter of 1-2 kpc, within which 
static thermal pressure dominates (Chevalier & Clegg 1985; Tomisaka & Ikeuchi 1986). The 
starburst region, i.e., the region producing the high pressure in this case, thus encompasses 
the two nuclei. The two nuclei are about 500 pc from the dynamical center of the system and 
probably have a velocity of about 100 km s _1 with respect to the dynamical center. Thus 
their crossing time across the central kpc is about 10 Myr. It is comparable to the typical 
timescale of a superwind. The two dominant energy sources thus move while a superwind 
develops, effectively making the region of energy injection larger than that of the two nuclei. 
Once a kiloparsec-size region of hot gas is formed by the moving energy sources, its evolution 
as a superwind is governed by the large scale structure of the ISM in which the hot gas is 
embedded. It is therefore conceivable that a main superwind centered around the dynamical 
center of the system blows out in the direction perpendicular to the main gas disk. This 
solution to the puzzle, however, seems to have a much room for improvement. 



6. Molecular Gas Properties 

6.1. temperature and filling factor 

Molecular gas in the center of NGC 3256 is warm. This can be seen in the map of peak 
brightness temperature shown in Fig. 4. The peak temperature in the map, in excess of 
10 K, is high considering that it is averaged over the 170 pc x 340 pc beam. The map is 
made by searching for the maximum brightness temperature in the velocity direction of the 
data cube at each sky position. The peak brightness temperature gives a lower limit on the 
temperature of molecular gas at each position, except for low intensity regions dominated 
by noise. It is a lower limit because of beam dilution — the filling factor of gas in the beam 
can be much smaller than unity — and also because the emission may be optically thin. In 
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particular, the map tends to give lower temperatures to regions with large line widths if the 
gas physical temperature and surface density are the same among regions. Despite these 
difficulties in reading the map, the map is indicative of higher gas temperature around the N 
nucleus and the tail region extending west from the S nucleus. As we noted, the N nucleus 
is the most active star forming region in the galactic center. The tail region corresponds 
to the tail feature seen in diffuse X ray emission (Lira et al. 2002, Fig. 12). The small 
temperature depression on top of the peak at the N nucleus and unremarkable temperature 
at the S nucleus are at least partly due to the dilution caused by the large line widths there. 

To put the beam diluted brightness temperature into perspective, we note that the peak 
value of 13 K in NGC 3256 (L 8 -iooo/im = 10 11 ' 56 L Q ; Sanders et al. 2003) is a factor of 3 
lower than that in the ultraluminous infrared galaxy Arp 220 (10 12 21 L Q ) observed with the 
same line at almost the same linear resolution (Sakamoto et al. 1999). It is also four times 
higher than the peak temperature of 3 K obtained in the starburst nucleus of M83 (10 1010 
L Q ) by convolving the CO (2-1) observations of Sakamoto et al. (2004) to match in linear 
resolution. All three galaxies are moderately inclined with respect to our line of sight (i ~ 
25° - 45°) and hence are expected to be similar in the degree of dilution due to galactic 
rotation. This comparison implies that molecular gas in galaxies with higher luminosity 
tends to be warmer or to have higher area filling factors or both when compared at 200 pc 
scale. 



6.2. gas cloud properties and line excitation 

6.2.1. line ratio observations 

NGC 3256 has a high 12 CO/ 13 CO intensity ratio. Figure 10 compares the spectra of 
the three CO lines obtained with a simulated 30" beam (FWHM) centered on the midpoint 
of the double nucleus. The flux densities integrated from 2590 km s _1 to 2970 km s _1 are 
4.62 x 10 3 , 1.37 x 10 2 , and 26.1 in unit of Jy km s -1 . The errors of these integrated intensities 
are ~10% due to flux calibration and ±2.6 Jy km s _1 due to noise directly estimated from the 
maps. The flux calibration error, however, cancels out in a line ratio from our simultaneous 
observations. Thus the error in the line ratio is dominated by noise. We obtain integrated 
temperature ratios of -Ri 2 /i 3 (2-l) = J( 12 CO 2- l)//( 13 CO 2- 1) = 30.9±0.6, i2i 2 /i 8 (2-l) = 
J( 12 CO 2 - 1)//(C 18 2 - 1) = 161 ± 15, and R 13/18 (2 - 1) = J( 13 CO 2 - 1)//(C 18 2 - 1) = 
5.2±0.5, where /(line) = fT h (lme)dv. Our R 12/13 (2-1) ratio agrees with those from previous 
single-dish observations (see Fig. 11). The it! 12 / 13 (2-l) ratio averaged over 5 kpc in NGC 
3256 is about three times higher than the ones observed in the central kpc of our Galaxy, 
10 ± 1 (Sawada et al. 2001), central 300 pc of the starburst galaxies IC 342 and NGC 253, 
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~ 9 (Meier et al. 2000, and Fig. 11), and in the central 500 pc of M82, ~12 (Mao et al. 2000). 
It is also higher than the -Ri2/i3(l~0) values of ~ 11 observed in the centers of nearby spiral 
galaxies (Sage & Isbell 1991; Aalto et al. 1991b; Paglione et al. 2001). On the other hand, 
the high i? 12 / 13 (2-l) ratio is within the range of values observed among luminous mergers, 
> 20 (e.g., Casoli et al. 1992b; Aalto et al. 1995; Glenn & Hunter 2001). 

The intensity ratios are not uniform in the galaxy. Our data indicate a trend that the 
i? 12 / 13 (2-l) ratio increases with the beam size (Fig. 11), suggesting higher ratios in the 
outer parts of the galaxy. The same trend is seen in i? 12 / 18 (2-l) with less significance. The 
i?i3/ig(2-l) ratio has too large uncertainty to detect its radial variation. The two nuclei 
have the lowest -Ri2/i3(2-l) value of 17 in our measurements. One can infer that the ratio 
is locally low at the two nuclei by looking at the integrated intensity maps in Fig. 2 where 
13 CO(2-l) shows peaks at the nuclei while 12 CO(2-l) does not. The radial gradient of the 
12 CO/ 13 CO intensity ratio in the J=2-l transition is opposite to that in the J=l-0 emission. 
The i?i2/i3(l _ 0) ratio at the center of the galaxy is high, 26-40 (Aalto et al. 1991a; Casoli 
et al. 1991; Becker & Freudling 1991; Garay et al. 1993), but it decreases to 11 at 43" NE of 
the galactic center (Aalto et al. 1995). 



6.2.2. gas model for the nucleus and the main disk 

The two-component model of molecular gas suggested by Aalto et al. (1995) seems to 
work for our higher central -Ri2/i3(2-l) value than in the centers of less luminous galaxies. 
In the model, molecular gas in a starburst nucleus has a significant amount of warm and 
turbulent envelope gas with moderate opacity (r[ 12 CO(l-0)] ~ 1 and r[ 13 CO(l-0)] <C 1) 
around cloud cores of high optical depths (r[ 12 CO(l-0)] ^> 1 and r[ 13 CO(l-0)] ~ 1). The 
large contribution (i.e., area filling factor) of the warm envelope gas makes the intensity- 
weighted optical depth of 12 CO emission moderate (~1) and makes the i? 12 / 13 (l-0) ratio 
higher than in core-dominated quiescent gas. The warm gas with high area-filling factor 
noted in §6.1 fits the profile of the envelope gas. The dense and high-opacity cores, on 
the other hand, are supported by the 13 CO(2-1)/ 13 CO(1-0) ratio close to unity, 1.3, in the 
center of NGC 3256 (Casoli et al. 1992b) and by the observations of dense-gas tracers such 
as HCN (Casoli et al. 1992a). For the .R12/13 ratio in the J=2-l transition, we note that 
the transition tends to have a higher opacity than the J=l-0 one. Thus we expect higher 
opacities in the envelope gas or, equivalently, a larger filling factor of optically thick cores. 
This and the presence of envelope gas make the i?i 2 /i3(2-l) ratio high at the galactic center 
but smaller than _R 12 / 13 (l-0), as observed. 

The model needs to be slightly modified in the outer disk to accommodate the opposite 
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radial trends in #12/13(1-0) and #12/13(2-1). For the outer disk, the same argument as 
above would lead to a lower #12/13(2-1) value in the disk than in the nucleus. It is because, 
according to the model, the disk emission is dominated by cloud cores. Such a radial trend 
of the J=2-l ratio contradicts our observations. One way to alleviate this is to assume 
that the CO excitation is subthermal in the disk cloud cores and that the cores have lower 
opacities in the J=2-l transition than they are expected to have from their J=l-0 opacities 
if the cores were in LTE. In other words, many cloud cores in the merger disk may look 
like low-opacity envelopes in the J=2-l lines. As a result, in the J=2-l lines, the bulk of 
emission in the galaxy has moderate (in 12 CO) and moderate to low (in isotopes) optical 
depths throughout, with opacities increasing toward the center. The radial gradient of the 
#12/18(2-1) ratio and the insignificant gradient of the #13/13(2-1) ratio are consistent with 
this. The subthermal excitation in the disk could be due to lower gas density there than in 
the nucleus, large velocity dispersion in the merged gas disk, and lower fractional abundance 
of CO with respect to H2. The higher [ 12 CO]/[ 13 CO] abundance ratio expected from the 
fresh supply of low- 13 C gas from peripheries of the merger progenitors (Casoli et al. 1992b; 
Henkel & Mauersberger 1993) would also increase the #12/13(2-1) ratio in low opacity clouds. 
Our conceptual model for the radial gradient of gas properties in NGC 3256 can be better 
constrained with more data. In particular, additional high-resolution data in other CO 
transitions, such as CO (1-0) or CO (3-2) would be most helpful to estimate the excitation 
conditions. 

Our observations in terms of radial gradient of the #12/13(2-1) ratio are consistent with 
those in the Antennae galaxy. Zhu et al. (2003) found that #12/13(2-1) was about twice as 
high in the overlap region of the merger than in the centers of the colliding galaxies. Their 
galactic center ratios, 13 ± 3 and 16 ± 4 for each galaxy, are marginally higher than the 
ones we quoted for local starburst nuclei. The higher #12/13(2-1) ratio in the overlap region, 
however, may be because the region has a starburst and resembles merger nuclei. Aalto 
et al. (1995) mentioned a decreasing #12/13(2-1) gradient from the nucleus outwards in the 
polar ring galaxy NGC 660 and an increasing gradient in the merging galaxy NGC 2146. 
Gas properties and excitation conditions of CO may be diverse among mergers in their outer 
disks. 



7. Discussion 

Our observations have provided a wealth of new information on the molecular interstellar 
medium in the center of NGC 3256. We can place them in the general context of the galaxy 
merger and that of merger- induced starburst. They also point areas for further studies. 
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The main gas disk is what is expected from the dissipative nature of the interstellar 
medium and what has been predicted from numerical simulations (e.g., Mihos & Hernquist 
1996; Barnes & Hernquist 1996). In the simulations, gas from the disks of progenitor galaxies 
collides and eventually coalesces into a single disk rotating around the dynamical center of 
the merger, with a strong concentration toward the center. Such a merged gas disk has been 
observed in a number of luminous and ultraluminous merging galaxies (e.g., Scoville et al. 
1997; Downes & Solomon 1998; Bryant & Scoville 1999; Tacconi et al. 1999). In line with 
those models and observations, the main gas disk of NGC 3256 has the majority of molecular 
gas in the system and the overall gas distribution peaks toward the center of the merger. 
The main gas disk is certainly supported by rotation, with the center of rotation located 
between the two nuclei. Presumably reflecting the violent past of the system, however, the 
disk shows signs of disturbances in the forms of short spiral arms and non-circular motions. 

Our observations reinforce the notion that properties of molecular gas in this merger as 
well as other luminous mergers are different from those in our Galaxy in the sense that the 
former contains more of warm, tenuous (i.e., high area- filling factor), and turbulent molecular 
gas (Aalto et al. 1995; Scoville et al. 1997; Downes & Solomon 1998). This must be due to 
the galaxy merging and to the starburst in the center of the merger. Our observations also 
suggest that gas properties change in spatial scales from over a kiloparsec to as small as 
a few 100 pc, which is our highest resolution. Examples of such spatial variation include 
the highest brightness temperature gas associated with the N nucleus that hosts the most 
active star formation in the galaxy, the local minima of 12 CO(2-l)/ 13 CO(2-l) ratio on each 
of the two nuclei, and the global radial gradient of the intensity ratio in the main gas disk. 
While these gas properties are probably the results of the merger and the starburst, the gas 
properties are also expected to influence the star formation from molecular gas. Therefore 
the study of gas properties at high resolutions helps to model the evolution of a starburst. 
We noted a trend of higher CO brightness temperature in more luminous galaxies. Statistical 
studies to probe trends like this and trends in line ratios will help characterize the molecular 
gas in mergers and starbursts. 

The molecular outflow from a merger has little observations in luminous or ultraluminous 
infrared galaxies of log(L IR / L e ) > 11.5, though surveys of atomic and ionized gas have shown 
that most luminous and ultraluminous infrared galaxies have a superwind (Heckman et al. 
1990, 2000; Martin 2005; Rupke et al. 2005). The recent detection by Takano et al. (2005) 
of high velocity gas in NH 3 absorption toward Arp 220 may be one of the rare observations 
of molecular outflow from a luminous infrared merger, although other interpretations are 
possible for the absorbing gas. In the model of superwind, an intense starburst deposits 
a large amount of kinetic energy from stellar winds and supernovae into the interstellar 
medium. The energy is thermalized in the ISM to generate hot gas of 10 6 -10 7 K that can 
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not be contained in the galactic potential and blows out of the galactic disk. Although 
molecular gas can not be hot enough to blow out by itself, it is expected to be entrained 
by the wind of hot gas. There are examples of such entrainment of molecular gas in nearby 
starburst galaxies of lower luminosities, e.g., in M82 (Nakai et al. 1987; Seaquist & Clark 
2001; Walter et al. 2002). It is therefore not surprising to detect a molecular outflow from 
such a luminous starburst as NGC 3256 in sensitive observations. The outflow parameters 
estimated in NGC 3256, i.e., the maximum outflow velocity, mass flux, momentum flux, and 
mechanical luminosity, are in order-of-magnitude agreement with those obtained from the 
survey of atomic Na absorption in luminous and ultraluminous infrared galaxies by Rupke 
et al. (2005). In particular, we reiterate the importance of gas outflow in the gas budget of 
the galactic center; molecular outflow can be nearly as significant as star formation in gas 
consumption. It is almost certain that we detected only a portion of the molecular outflow 
in NGC 3256, limited by our sensitivity. Observations at higher sensitivity and resolution 
would shed more light on the mechanism for molecular gas entrainment and the effects of 
the evacuation of molecular gas on the starburst evolution. Such observations will also help 
testing the alternative interpretations of the high velocity gas. 

Regarding the double nucleus, we suggest that the gas peaks around and the steep 
velocity gradient across each nucleus are most likely due to a mini disk rotating around the 
nucleus, although the mini disks certainly merit further observations for confirmation. There 
is a striking similarity among NGC 3256, Arp 220, and NGC 6240 in terms of their twin 
(possible) remnant nuclei of a few 100 pc sizes and ~10 9 M Q masses within the projected 
separation of 1 kpc. In the earlier phase of galaxy collision, numerical simulations suggest 
that a sizable fraction of gas in each galaxy can be funneled to the nucleus of each galaxy, 
forming a massive rotating core of gas and that of stars through starburst (Noguchi 1988; 
Barnes & Hernquist 1996; Mihos & Hernquist 1996). Observations of interacting galaxies 
find such starbursting gas concentrations in many, though not all, galaxies (Aalto et al. 
1997; Wilson et al. 2000; Yun & Hibbard 2001; Kaufman et al. 2002; Evans et al. 2002). As 
the merger proceeds, stellar systems will merge from outside, leaving nuclear stellar cores 
until the final stages (Funato et al. 1992; Boylan-Kolchin & Ma 2004). Gas from the merger 
progenitors will also coalesce first from outside to form a large merged disk, which is most 
likely the main gas disk in NGC 3256. In the mean time, the nuclear gas trapped deep inside 
the gravitational potential of each progenitor likely remains within the nucleus. The pairs 
of massive nuclei in the three merging galaxies are plausibly the remnants of such massive 
cores of stars and gas. (Gas was presumably striped from the stellar cores in NGC 6240.) 
The massive remnant nuclei should strongly perturb the larger, merged gas disk as they sink 
toward the dynamical center of the system owing to dynamical friction. Recently Escala et 
al. (2005) suggested from numerical simulations of a binary massive black hole in a merger 
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that cold gas in the central merged disk can be a significant source of dynamical friction and 
that the gas disk is affected in return forming a gap or spirals or both around the binary. 
The situation may be similar in NGC 3256. The gas concentration around each nucleus 
should eventually contact the gas around the other nucleus or in the main disk, loses angular 
momentum, and probably falls toward the dynamical center of the system to become a part 
of the main gas disk. The complicated dynamics of gas, stars, and presumably massive black 
holes, remains to be studied along with the evolution of starburst in the perturbed region. 
Our observations show that NGC 3256, along with Arp 220, is one of the best targets for 
such studies. 



8. Conclusions 

We have made the first interferometric observations of molecular gas and dust in the 
center of the southern luminous merger NGC 3256, improving the spatial resolution over 
previous single-dish observations by an order of magnitude. Our main findings are the 
following: 

1. There is a large disk of molecular gas (r > 3 kpc) with a strong concentration toward 
the double nucleus. The velocity field of the disk suggests rotation around a point 
between the two nuclei. The mass of neutral gas in the central 20" (r < 1.7 kpc) is 
estimated to be ~3 x 1O 9 M . Our observations suggest that the molecular gas in the 
merger system has mostly merged and settled into this main disk at the radii of 1-5 
kpc. 

2. In the central kiloparsec, our highest resolution (1" x 2" = 170 pc x 340 pc) image 
shows that molecular gas peaks near or at each of the two nuclei. There is also a steep 
velocity gradient across each nucleus, suggestive of a small gas disk rotating around 
the nucleus. The dynamical mass of each nucleus is estimated from the gas velocity 
to be ~10 9 M Q within a radius of 1" (= 170 pc). The suggested configuration — two 
mini disks around the binary nuclei along with a larger main disk — is similar to that 
in Arp 220. 

3. We discovered high velocity molecular gas of ~10 7 M Q at the galactic center between 
the two nuclei. It is up to about 400 km s _1 from the systemic velocity of the galaxy. 
The terminal velocity is about two times larger than that due to rotation of the main 
gas disk. We have modeled the high velocity gas as a molecular outflow, entrained by 
the superwind that the starburst in the central kiloparsec is known to have created. 
Outflow parameters, including the outflow rate of ~10 M Q yr -1 , are estimated from 
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the data. They suggest that the dispersal of molecular gas from the starburst region 
can be as significant as the gas consumption by star formation in the gas budget. 

4. Our simultaneous observations of the 12 CO(2-l), 13 CO(2-l), and C 18 0(2-l) emission 
provided accurate measurements of intensity ratio among those lines. The 12 CO/ 13 CO 
intensity ratio in the J— 2-1 transition is higher at the nucleus than in the centers 
of our Galaxy and nearby starburst galaxies. Our observations are consistent with 
the model that the intense starburst in the merger nucleus is taking place in warm 
turbulent gas with a high filling factor and moderate 12 CO opacity. Positive radial 
gradients of 12 CO/ 13 CO and 12 CO/C 18 ratios are observed in the J— 2-1 transitions. 
They may reflect subthermal excitation in the merger disk. 
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Fig. 1. — Optical and near-IR images of NGC 3256. (left) Blue-band image from the UK 
Schmidt Survey through the Digitized Sky Survey. The image is 10' by 6' (= 100 kpc x 60 
kpc). (right) i^-band image of NGC 3256 taken by Rothberg & Joseph (2004). The central 
2' by 2' (= 20 kpc x 20 kpc) is shown. In each image, north is up and east is to the left, 
and intensity is shown with logarithmic scale. 
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Fig. 2. — Natural-weighting maps of CO line and continuum emission in the center of NGC 3256. 
(a) 12 CO(2-l) integrated intensity. The n-th contour is at 7.1 x n L75 Jy beam -1 km s -1 . (b) 
12 CO(2-l) mean velocity (i.e., 1st moment). The velocity contours are in steps of 25 km s _1 . (c) 
13 CO(2-l) integrated intensity. Contours are at [1,3,5,7,9] x 2.1 Jy beam" 1 km s" 1 . (d) 13 CO(2- 
1) mean velocity, (e) C 18 0(2-l) total intensity map, integrating from 2650 km s" 1 to 2910 km s" 1 . 
Contours are in steps of 3.2 Jy beam -1 km s _1 (= 2a). (f) 1.3 mm (222.7 GHz) continuum map 
made by combining the data from both sidebands. Contours are in steps of 5.0 mJy beam" 1 (= 2a). 
In each panel, the double nucleus is marked with a pair of crosses. The one in the north is the 
nucleus N and the other in the south is the nucleus S. The synthesized beam of each map is shown 
at the bottom left corner. 
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Fig. 3. — Inclination-corrected radial distribution of 12 CO(2-l) emission. The mean inte- 
grated intensities (left axis) are measured in concentric rings in the galaxy plane, for which 
we adopt an inclination of 45°, a position angle of the major axis of 90°, and the center at the 
midpoint of the two nuclei. The geometry is consistent with the kinematical fit of our CO 
velocity data. The cumulative flux (right axis) is calculated from the integrated intensities. 
The line flux is measured from the integrated intensity map made by summing up channel 
maps, rather than a moment map. The natural weighting dataset is used for this plot after 
correcting for the primary beam attenuation. The flux scale has a ~10% uncertainty. 
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Fig. 4. — High-resolution maps of 12 CO(2-l) emission in the center of NGC 3256. (a) 
Integrated intensity in unit of Jy beam -1 km s™ 1 . (b) Mean velocity in unit of km s -1 . (c) 
Line width (i.e., the 2nd moment) in unit of km s . (d) Peak brightness temperature in 
unit of K. The two nuclei of the merger are marked with + signs. The synthesized beam is 
at the bottom left corner. The n-th contour in the integrated intensity map is at 3.9 x n 1 ' 5 
Jy beam -1 km s _1 . The velocity contours are in steps of 25 km s -1 and 10 km s _1 in the 
1st and 2nd moment maps, respectively. The temperature contours are in steps of 1.5 K 
starting from 3 K. The temperature map, only to which the primary beam correction is 
applied, is masked in its noisier outer region. The high-resolution data cube was made with 
super-uniform weighting and recovered 54 % of the total CO flux in our natural weighting 
data shown in Fig. 2. 
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Fig. 5. — SMA CO (2-1) contours overlaid on an HST /-band image of the center of NGC 
3256. The CO contours are the same as in Fig. 4, and the optical image is shown with 
logarithmic scale. The radio positions of the two nuclei are marked with plus signs. The x 
between the two nuclei is the dynamical center of the galaxy estimated from the CO velocity 
field. North is up, and east is to the left. This figure covers an area of 34" by 34". See Zepf 
et al. (1999) for the HST observations and a color composite image of the same area. 
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Fig. 6. — Position-velocity (PV) diagrams made from the high-resolution data cube, (a) PV 
diagram along P.A.= 270° across the N (north) nucleus, (b) PV diagram along P.A.= 270° 
across the S (south) nucleus, (c) PV diagram along P.A.= 270° across the middle point (M) 
of the double nucleus, (d) PV diagram along P.A.= 0° across the two nuclei. Each PV cut 
has a 2" width except the one across M that has a Of! 2 width to minimize the contamination 
by the gas around the N and S nuclei. The black rectangle at the bottom-left corner of each 
panel shows the spatial resolution along the PV cut and the 10 km s _1 velocity resolution of 
the data. Contours are in steps of 97 mJy beam" 1 or 1.2 K in the brightness temperature. 
The dotted line in each diagram shows the systemic velocity of the galaxy (V^ ys = 2775 
km s _1 ). 
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Fig. 7. — Position-velocity diagram along P.A.= —90°. The dotted line indicates the systemic 
velocity of the galaxy for which we adopt V r sys (radio, LSR) = 2775 km s _1 . The zero point 
of the horizontal axis is where the double nucleus is. Intensity is in Jy beam -1 and can be 
converted to brightness temperature in kelvin by multiplying 2.8. The ra-th contour is at 
n 1 ' 75 x 30 mJy beam -1 . The black rectangle at the bottom left corner tells the spatial and 
velocity resolutions. The arrows point the high velocity gas in the galactic center. 
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Fig. 8. — Molecular gas at extreme velocities. The maps average velocities where either 
redshifted (a) or blueshifted (b) high velocity gas is seen in the galactic center. Contours 
are in steps of 2 a in each map, i.e., 16.6 mJy beam -1 for the redshifted map and 24.6 mJy 
beam -1 for the blueshifted map. The plus signs mark the two radio nuclei. The synthesized 
beam is on the bottom left corner of each panel. The arc-line emission 15" west of the nucleus 
in the blueshifted map is not the high velocity gas but corresponds to a spiral feature in the 
outer disk, where noncircular motion is significant. 
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Fig. 9. — Spectrum of 12 CO (2-1) emission in the central 5" of NGC 3256. The bottom frame 
is to show the wing emission from the high velocity gas. The spectrum is the one at the 
midpoint of the binary nuclei. The r.m.s. noise in this profile is 30 mJy beam -1 ; the noise 
was measured in the image domain in the data cube from which the spectrum was obtained. 
The spectrum is from the natural- weighing data cube convolved to a 5" resolution (FWHM). 
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Fig. 10.— Spectra of the three CO lines in the central 30" of NGC 3256. The natural- 
weighting data cube of each line was corrected for primary-beam attenuation and then con- 
volved to the 30" resolution (FHWM). The spectrum was then measured at the position of 
the phase tracking center, which is the middle point of the double nucleus. The dashed line 
in the 12 CO panel is the same 12 CO spectrum multiplied by 10. 
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Fig. 11. — Intensity ratio of the J =2-1 transitions of CO isotopomers in the centers of 
NGC 3256 (open symbols) and NGC 253 (filled symbols). Circles, squares, and diamonds 
are for 12 CO/ 13 CO, 12 CO/C 18 0, and 13 CO/C 18 0, respectively. The ratios in NGC 3256 
are computed between brightness temperatures integrated over the range of 2590 km s _1 - 
2970 km s -1 . The rightmost and middle data of NGC 3256 are measured at the midpoint 
of the double nuclei from the natural weighting data convolved to 30" and 10" resolutions 
(FWHM), respectively. The left group of NGC 3256 data are measured at the N nucleus (left 
points) and S nucleus (right points) in the natural weighting data convolved to 5"2 x 2"0 
resolution. The nuclear data are assigned the beam size of 3.2 (=\/5.2 x 2.0) arcsec and are 
given a slight horizontal offset for each nucleus to avoid overlap. Correction for the primary 
beam attenuation was applied before each convolution. Error bars are for ±la. The two 
12 CO/ 13 CO data points for the beam size of about 25" are from the single dish observations 
of Casoli et al. (1992b) and Garay et al. (1993). The data for the center of NGC 253 are 
from the SMA observations of Sakamoto et al. (2006) convolved to a 250 pc = 15" (FWHM) 
resolution and integrated over the full velocity range of 50 - 450 km s -1 . (The beam size of 
250 pc is l'/5 at the distance of NGC 3256.) The la errors of the NGC 253 data are smaller 
than the symbols. 



